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Études structurales en forêt 
par les techniques de morphologie mathématique 


J. BACHACOU *. J. BOUCHON ** et S. TOMIMURA ** 


RÉSUMÉ 


Cette étude rassemble et prolonge des travaux présentés par ailleurs (BACHACOU. 1978: BOUCHON, 1977; 
BOUCHON, 1975). Elle présente des techniques permettant de quantifier les impressions visuelles reçues par 
l'œil lorsqu'il observe des ensembles de taches : projections horizontales des couronnes d'un peuplement de 
Sapin ou carte des taches d'une régénération naturelle. 


SUMMARY 


This study is a prolongation of precedent works (BACHACOU, 1978; BOUCHON. 1977; BOUCHON, 197$). 
The techniques presented here are based on the method of regionalized variables and mathematical 
morphology (SERRA. 1969: MATHERON, 1970). The initial images are simplified towards images or indexes 
easier to analyze. Applications are presented concerning the map of horizontal projection of fir crowns 
(Abies alba) and the map of a natural regeneration of beech (Fagus silvatica). 


Les méthodes présentées ici visent principalement à familiariser les lecteurs avec 
des procédés encore peu connus. En ce qui concerne la théorie, nous renvoyons à 
J. Serra (1969): pour la technique plus particulière d'analyse des images: il est conseillé 
de voir J. Serra (1967) et BacHacou (1978). 


1. SYMBOLES EMPLOYÉS 


Ri; Rio. point courant de la grille appliquée sur la carte: 

T. taux de régénération: 

M. M. nombres utilisés pour le calcul de la précision des estimations de superficies par comptage de 
points: 

Max, maximum (AM, Mt 

min. minimum (M,, Mt 

M. nombre total de points R; tombant dans la régénération: 

L, longueur totale des transects; 

H. somme des longueurs des segments de transects intérieurs à la régénération: 

N, nombre de taches traversċes par les transects; 
(Li. Le. Hi, He, Ni, No. valeurs correspondantes pour les lignes et les colonnes): 

N (u), nombre de segments de longueur u rencontrés par le transect; 

€. périmètre spécifique; 

H. traversée moyenne=H/N =(H, HUN, +N): 

h, diamètre moyen des taches =4 H/x:; 


* C.NR.F./I.N.R.A. Station de Biomėtrie, Champenoux. 54280 Seichamps. 
** C.N.R.F./LN.R.A., Station de Sylviculture et de Production. Champenoux, 54280 Seichamps. 
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surface moyenne des taches: 

nombre de segments par unité de longueur: 
nombre de taches par unité de surface: 
granulométries à une et deux dimensions: 
longueur des segments: 

diamètre des taches: 

érosions et dilatations; 

pourcentage de recouvrement de deux strates: 
rayons des étoiles 


2. PEUPLEMENTS FORESTIERS ÉTUDIÉS 


L'un des peuplements étudiés est décrit dans BoucHon (1973) : il s'agit d'une 
réserve biologique intégrale. dite La Tillaie, en forêt domaniale de Fontainebleau, en 
parcelles 270 et 271: ce peuplement est situé sur un plateau entre 133 et 139 m 


N REGENERATION 


FIG. 1. — Forêt Domaniale de Fontainebleau : peuplement de « La Tillaie » 


Carte de la régénération. jusqu'aux stades du fourré et du gaulis fermés. 


FiG. 2. — Forêt Domaniale de Fontainebleau : peuplement de « La Tillaie ». 
Carte des sols : en noir. podzols et sols podzoliques: en gris : sols néopodzoliques. 


Œcologia Plantarum 


ÉTUDES STRUCTURALES EN FORÊT 207 


Bas et Hauts 
PERCHIS 


FiG. 3. — Forêt Domaniale de Fontainebleau ` peuplement de « La Tillaie » 
Carte des bas et hauts-perchis. 


d'altitude. La température moyenne annuelle est de $.85°C et la pluviosité de 696 mm. 
Le substrat est formé en surface par le « limon des plateaux » reposant sur le calcaire de 
Beauce en majorité. et en partie au nord-ouest sur le grès et les sables de Fontainebleau. 
Le peuplement forestier est presque exclusivement constitué par une hêtraie contenant 


Movenne et 


Haute FUTAIE 


FIG. 4. — Forêt Domaniale de Fontainebleau : peuplement de « La Tillaie ». 
Carte de la moyenne et haute futaie. 


quelques vieux Chênes disséminés qui ne se régénèrent pas; les taches de régénération 
sont donc uniquement constituées de jeunes Hêtres : voir les figures 1 à 4 pour les 
cartes de la régénération, des sols, des bas et hauts perchis et de la moyenne et haute 
futaie. 
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Pio 5. — Projection horizontale des couronnes dans un peuplement naturel de Sapins : Abies alba Mill. 
(B. Vins. 1961). La granulométrie à une dimension est obtenue en mesurant les longueurs des cordes 
interceptées par les projections des couronnes le long de divers transects. Ici on a fait sept transects dans 
une direction (1 à 7) et onze dans l'autre (4 à K). Les mesures faites le long du transect n° 3 ont été 
figurées. 


Le deuxième peuplement est une sapinière décrite dans Vins (1961). Il s’agit d’un 
peuplement naturel d'un demi hectare, situé à 1025 m d'altitude, dans la Haute Tatra; 
la pluviosité annuelle est d'environ 1 000 mm et la température moyenne de 4°C. Les 
données consistent en la projection horizontale des cimes (voir fig. 5). 


3. MESURES DE SUPERFICIES PAR COMPTAGE DE POINTS 
… Le but était d'étudier la régénération dans la réserve naturelle de La 
Tillaie (fig. 1). Une grille de points a été appliquée sur la carte et un simple comptage 
indique que la surface régénérée représente 1=28.28 °, de la superfcie totale. Une 
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méthode indiquée dans Bouchon (1975) permet de calculer la précision de cette 
estimation; la grille appliquée sur la carte définit des lignes et des colonnes de points; il 
est nécessaire de calculer deux nombres : 

— M, pour les lignes; 

— M, pour les colonnes. 


Règle pour le calcul de M. 


Lorsqu'on passe de la colonne i à la colonne i+ 1, on peut définir quatre cas pour 
deux points R; et R;,, adjacents : 

a) R;et R;,, appartiennent tous les deux à la régénération; 

b) R;et Ri+ı n'appartiennent pas à la régénération; 

c) R; appartient à la régénération, mais R;., ne lui appartient pas: 

d) R; n'appartient pas à la régénération, mais R;,, lui appartient. 

La somme des cas c et d constitue la « variation diamétrale » entre la colonne iet la 
colonne i+ 1. La somme de toutes ces variations diamétrales intercolonnes obtenue en 
balayant le champ est égale à M.. On obtient M, de la même manière pour les lignes. 

Si Max est la plus grande des deux valeurs M, et M. et si min est la plus petite, 
alors la variance de l'estimation de T est donnée par (MATHERON, 1970) : 


M étant le nombre total de points tombant dans la strate régénérée. Les calculs 
montrent qu'au seuil de 0,05 la surface régénérée représente 25,5 à 31,1 % de la 
superficie totale. 


4. MESURES DE SUPERFICIES PAR TRANSECTS 


Pour faire cette mesure, il suffit de tracer des lignes balayänt le champ 
d'observation et de calculer le rapport entre la somme H des segments intérieurs aux 
taches de régénération et la longueur L totale du transect. 


On a pris comme transects les droites joignant les points ayant servi à faire 
l'estimation du paragraphe précédent. 
Transects parallèles aux lignes 


La somme des traversées intérieures aux taches est H, =2293 mm (sur la carte). 
La somme des transects est L , =7920 mm. D'où un taux de régénération de : 


2293/7920= 28,95 %. 


Transects parallèles aux colonnes 
Taux de régénération : 


H./L.=2339/7920=29,53 %. 


Ces deux valeurs sont compatibles avec le résultat du paragraphe précédent. Dans 
la suite, nous prendrons : 


T=(H;+H.)/(L;+L.)=2924% 


comme estimation du taux de régénération. 
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5. RETOUR AUX PROBLÈMES DU FORESTIER 


Le simple taux de régénération n'est pas un critère suffisant pour le forestier; il est 
nécessaire de connaître la répartition spatiale des taches, leur forme, leur taille 
unitaire... pour expliquer les mécanismes de la régénération, pour décider de l'utilité 
éventuelle d'opérations sylvicoles visant à compléter le peuplement, etc. Les 
paragraphes suivants présentent chacun un essai d'analyse morphologique appliquée 
aux ensembles de taches qui sont représentées par les figures 1 à 5. 


6. RÉGULARISATION DES IMAGES 
Un appareil, l'analyse de texture (SERRA, 1967; BacHacou. 1978) permet de 


simuler une opération qui intéresse vivement le forestier. Si on observe attentivement la 
figure 1, on y voit de toutes petites taches de régénération qui ont peu d'avenir si elles 


f 7 


FIG. 6. — Le ouverture » est une opération qui consiste à éliminer d'une image les petites taches. La 
« fermeture » est l'opération inverse ` elle élimine les petites lacunes: ces opérations ont été appliquées à 
la régénération de la figure 1: on a successivement fait une fermeture d'amplitude 10 m (image du haut) 
qui a éliminé les trouées inférieures à 10 m, puis une ouverture de mème amplitude (image du bas). Noter 
que ces deux opérations ne sont pas commutatives. 


restent isolées, car la morphologie des arbres croissant librement est mauvaise du point 
de vue sylvicole: inversement. les petites lacunes intérieures à de grandes taches de 
régénération se refermeront au fur et à mesure que le peuplement grandira: l'analyseur 
de texture, par des opérations successives dites fermeture et ouverture, permet de 
gommer automatiquement les petites trouées et les petites taches: la figure 6 ainsi 
obtenue est plus conforme à la réalité sylvicole, en ce sens qu'elle montre clairement au 
forestier où il lui reste à porter son effort. 
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7. SUPERPOSITION DE DEUX CARTES 


Les deux premières cartes (régénération et sol : fig. Let 2) montrent une relation 
évidente entre la présence des podzols et des sols podzoliques et l'absence de 
régénération. Devant la banalité du résultat qu'aurait donné l'étude simultanée de ces 
deux cartes, l'étude de costructure n'a pas été faite ici UL 

MaRBEAU (1973) a voulu comparer les cartes des bas et hauts perchis d'une part et 
de la moyenne et haute futaie d'autre part. Pour cela, il a érodé puis dilaté les taches de 
futaie et calculé les teneurs en perchis dans les couronnes successives parallèles aux 
limites de la futaie. En appelant A (— i) les érosions et A ( + i) les dilatations, les teneurs 
obtenues sont les pourcentages de perchis (par rapport à la superficie de chaque 
couronne). On obtient les résultats indiqués dans le tableau I. 


TABLEAU I 


Pourcentage de perchis au voisinage des taches de futaie: on a calculé ces teneurs dans des bandes 
de 4m parallèles aux limites de la furaie, à l'intérieur pour les indices négatifs, à l'extérieur 
pour les indices positifs 


La teneur moyenne du perchis sur l'ensemble de la parcelle est de 25 ",. On peut 
conclure à une certaine indépendance entre la présence de la futaie et celle des perchis 
[l'érosion A(—4) amenant presque au cœur des plus grosses taches de futaie]; cette 
méthode conduirait à affiner le concept d'essence d'ombre /essence de lumière (°): une 
étude analogue serait intéressante dans un peuplement voisin où la proportion de 
Chêne (essence de lumière) est importante; la quasi-indifférence entre P;et la distance à 
la lisière des futaies conduit à classer le Hêtre parmi les essences d'ombre. La faible 
valeur de P,=12 °,, peut s'expliquer par le fait que les faines étant des graines lourdes. 
en l'absence de semenciers. la présence de jeunes peuplements est impossible. 


8. PÉRIMÈTRE SPÉCIFIQUE 


C'est, ramenée à l'unité de surface, la longueur de la ligne de contact entre deux 
strates. On peut avoir très facilement une estimation du périmètre spécifique dans le cas 
des taches de régénération de la figure l. à l’aide des transects déjà utilisés au 
paragraphe 3. En effet. si L est la somme des longueurs des transects et si N est le 
nombre de taches traversées, on montre que le périmètre spécifique peut être calculé 
par la formule : 


C=rN/L. 
Ici, on a pris : 
N=Ni+N. et L=L;+L,, 


(Ur En réalite, la superposition avec la carte de la végétation qui n’a pas ète reproduite ıcı était encore 
plus explicite : la régénération est absente des taches de fougère aigle. 


(È) Chez les essences d'ombre. la régénération et la croissance à « l'ombre » du peuplement dominant 
sont possibles. Chez les essences de lumière, le maintien du couvert empèche l'installation de la régénération. 
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sommes des valeurs obtenues suivant les lignes et les colonnes et on a trouvé : 
C=0,0506m" ", 


ce qui signifie que dans un hectare, on a 506 m de contact entre les taches de 
régénération et les vides. 


Ce périmètre permet de calculer un indice de dispersion des taches. En effet, en 
admettant que les taches de régénération sont circulaires et égales on peut calculer 
combien il faudrait de taches pour qu'elles couvrent la superficie T (calculée au 
paragraphe 2) en ayant comme périmètre spécifique total C. Un calcul simple montre 
que le nombre de taches circulaires qu'il faudrait est donné par : 


8=C?/(4rT), 


8 étant un nombre par unité de surface. 
Avec C=0,0506m"! et T=0,292 4 on trouve : 


8—6,927 taches/ha, 


ce qui, étant donné le taux de régénération de 0,2924 conduit à une dimension 
moyenne des taches de S —420 m°. 


En réalité, l'observation de la figure 1 montre que la régénération n'est pas 
distribuée uniformément et que les taches n'ont pas une forme circulaire: les deux 
valeurs : 


8—6.97 taches/ha et  S—420m° 


ne sont donc que des indices de dispersion: ils n’ont aucune réalité biologique. 


9. GRANULOMÉËTRIE DES TACHES (exemple : cimes de sapins) 


Supposons que sur la figure 5 on mesure les diamètres des couronnes qu'on 
assimile à des cercles. L'ensemble de ces diamètres présente une certaine distribution: 
on peut définir la densité de couronnes de diamètres supérieurs à h par la fonction : 


8:(1—F(h) 
où 6, est le nombre moyen de couronnes par unité de surface et où F;(h) est appelé 


granulométrie à deux dimensions. 


De la même manière, un transect tel que celui qui est représenté sur la figure 5 
permet de calculer une distribution des traversées: on définit la densité des longueurs 
traversées supérieures à u par la fonction : 


6, (1—F,(u)) 
où 0, est le nombre moyen de segments par unité de longueur et où F, (u) est appelé 


granulométrie à une dimension. 


On montre (MATHERON, 1970) que si les taches sont circulaires, ou voisines de 
cercles, on peut calculer la distribution à deux dimensions par la simple connaissance 
de la distribution à une dimension : 


28 
82(1— F2 (h)) af 


= Filu) 


flah? 
t AHB 


du. 
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Dans une étude sur les couronnes de Sapin de la figure 5 (BOUCHON, TOMIMURA, 
1977), on a montré une très bonne coïncidence entre la distribution 8,(1 — F, (h)) 
mesurée directement sur la carte et celle calculée à partir de cette formule (*). Mais la 
même méthode appliquée aux taches de régénération ne peut rien donner, car ces 
taches ne sont pas circulaires, ni même convexes; cependant, des indices analogues à 
ceux calculés au paragraphe précédent peuvent être calculés à partir de la 
granulométrie à une dimension. 

Par exemple, on peut remarquer que la connaissance de la granulométrie à une 
dimension permet un calcul direct de la valeur S correspondant à la surface moyenne 
des taches. On peut montrer que si on balaye un cercle de diamètre h par des lignes 
parallèles, la valeur moyenne des longueurs des cordes interceptées est : 


H=rh/4. 


Inversement, la granulométrie à une dimension permet le calcul direct de la 
traversée moyenne: ici, on aura par exemple ` 


H=(H;+H)/N;+N)=18.16m. 


Le diamètre moyen h correspondant serait 23.1 m, ce qui correspond bien à un 
cercle de surface 420 m°. Mais soulignons encore que ces valeurs ne sont que des 
indices. puisqu'elles reposent toutes sur l'hypothèse que les taches sont circulaires et de 
mème taille. 


Un meilleur indice peut être calculé à partir de l'intégrale : il ne fait plus 
l'hypothèse que les taches sont toutes égales. 


On peut en effet calculer : 
20, | ° Fi(u) 
m Jo u 


du. 


CE 


ce qui devient, en partant directement de la distribution des longueurs des segments : 


2 Z Nu 


L étant la longueur totale des transects (ici L= L ,+L,.=2 x 7920 mm), N (u) étant le 
nombre de segments de longueurs u. 
On trouve H. 12.58 taches de régénération ha. 


Enfin, en ajustant une courbe 8, (1 — F, (u)) sur la distribution des longueurs des 
segments interceptés par les transects et en calculant l'intégrale 0, (1— F (h)), on peut 


avoir une estimation de la surface moyenne des taches de régénération S et du taux de 
régénération T. 
On trouve : 


S=185m?, T—2323%. 


D'où les valeurs indiquées dans le tableau récapitulatif pour h, Cet H 
(voir tableau II). 


(°) Voir en annexe pour l'utilisation pratique de cette formule. 
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TABLEAU I 


Comparaison des valeurs trouvées par les diverses méthodes. La plupart de ces valeurs peuvent être 
calculées sans avoir à formuler d'hypothèses sur la forme des taches: cependant : les résultats 
suivis d'une astérique supposent que les taches sont convexes: les résultats suivis de deux 
astériques Supposent en plus que les taches sont de même taille. 


Par comptage 


Par transect Granulométrie 
de points 


Traversée moyenne : H 18.20 m KA 


Nombre moyen de taches par | 
unité de surface : 8, 6.97 taches/ha** | 12.58 taches/ha* 


Périmètre spécifique : C 506 m/ha 605 m/ha* 


Diamètre moyen des taches : 23. mes 15.3 m* 


Surface moyenne des taches : 420 m2** 


Taux de régenération : T...... 28. 39.24 °, 
ech 


10. ÉTOILE DES TACHES 


On a vu que dans le cas de taches non circulaires, ou plus généralement non 
convexes, les indices qu'on calculait par les méthodes de granulométrie sont difficiles à 
interpréter. I] existe un paramètre de forme adaptée aux formes non convexes : c'est 
l'étoile. 

Définition : supposons une image telle que la figure 7. sur laquelle on peut définir 
deux strates : une strate de taches et une strate de vides. Un point quelconque R se 
trouve à l'intérieur d'une des deux strates: l'étoile de R est le lieu des points 3 tels que 
les segments RY soient entièrement contenus dans la strate à laquelle appartient R. 

On voit que pour une tache convexe, quel que soit le point R, intérieur à cette 
tache, l'étoile de R, est identique au grain tout entier. Unien est plus ainsi quand les 
taches ne sont pas convexes. 


Estimation de l'étoile moyenne des taches 


La définition de l'étoile d'un point montre que si l'on connaissait l'équation 
polaire r(w) de son contour relativement au point R pris comme origine des 
coordonnées, alors on pourrait calculer la surface de l'étoile E par l'équation : 


E=- j r? (w) do. 
J 0 


Si l'on considère le point R; et si l'on mesure d rayons rij de son étoile faisant entre 
eux des angles de 27/q. on montre que ` 


est une estimation non biaisée de E. 
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On a ainsi pour l'ensemble des points R; intérieurs aux taches de régénération 
mesuré huit rayons décalés de x/4 et estimé les surfaces Ê; des étoiles des R;. 


On trouve : 
— ZE 157070 
= TT = —34]5 mi. 
E M 46 3415 m 


Cette surface est très grande par rapport aux estimations précédentes: mais, étant 
donné le mode de calcul qui privilégie les grandes étoiles, ces estimations ne sont pas 
réellement comparables: en effet, le nombre de points échantillons R; tombant dans 


Fi, 7. — Sur une image composée de plusieurs strates (ici on observe une strate de grains et une strate de 
vides), on définit l'étoile d'un point R comme l'ensemble des points visibles de R et appartenant à ia 
même strate que R. 


une tache est proportionnel à la surface de cette tache: pour avoir une estimation de la 
surface moyenne des taches par la méthode des étoiles. il faudrait donc pondérer 
chaque mesure précisément par son inverse : 


_DU/ÉJÉ AM 


R — =354 m°. 
Fia "EBI D 


La comparaison de cette valeur à celles obtenues précédemment pour la surface 
moyenne des taches de régénération est instructive : 

a) si les taches étaient convexes, la superficie moyenne calculée à l’aide des étoiles 
serait égale à la surface moyenne des taches calculée par la granulométrie. La différence 
entre les deux valeurs (354 et 185 m°) permet de chiffrer l'écart à la convexité: 

b) compte tenu du taux de régénération. le périmètre spécifique de 506 m/ha 
indique que la forme des taches de régénération est très irrégulière: pour obtenir un tel 
périmètre, il faut : 

— soit beaucoup de petites taches, 

— soit des taches présentant de nombreuses ramifications et des isthmes étroits: 

c) or. la traversée moyenne et la surface moyenne indiquées par les mesures de 
transects sont importantes. Il faut donc en déduire que la valeur trouvée pour le 
périmètre spécifique correspond plutôt à des taches très ramiñées qu'à un grand 
nombre de petites taches: cela confirme bien l'impression qu'on a en observant la 
figure 1. 
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CONCLUSION 


La caractérisation des formes est un problème complexe. Le plus souvent seule la 
comparaison de plusieurs des indices présentés ci-dessus permet de décrire 
synthétiquement les structures spatiales observées. 

En toute rigueur : 

— la méthode des transects ($ 4 et 8) n'est utilisable que lorsque les taches sont 
circulaires et de même taille: 

— l'étude des granulométries ($ 9) n'est possible que si les taches sont circulaires: 

— l'estimation de la surface moyenne par les étoiles n'est possible que lorsque les 
taches sont convexes. 

Dans le cas général, que les taches soient non convexes ou qu'on s'intéresse aux 
vides au lieu des taches, l'étude des étoiles est la seule méthode ne présupposant rien sur 
l'allure des formes. 


ANNEXE 


CALCUL DE LA GRANULOMÉTRIE A DEUX DIMENSIONS 
A PARTIR DE LA GRANULOMÉTRIE A UNE DIMENSION 


La formule qui permet de passer de l’une à l'autre est : 


Flu) 


Se, 


db H ft 


= 20, 
DÉI F: th= — 
R 


où 8;(1 — F; (h) est la densité de grains de grosseurs supérieures à h: H. est le nombre moyen de grains par unitè 
de volume (surface ou longueur): F;(h) est la granulomètrie à į dimensions. 


1. On détermine expérimentalement la courbe 8, (1—F, ut, 

2. On ajuste cette courbe à un polynôme de degré supérieur ou égal à 2 pour que F (0) soit définie. 
3. On calcule 6, F; (u). 

4. On calcule l'intégrale en remarquant que ` 


pp page SE VS een 


ur 2 pi 
A = La ees 
`Hh p+] ka pip—2)...(p—2i) 


avec 
FE@m+1)=0. 
sip=2m+l LE 
l —— =m, 


2n—1)(2n—3)...1 ES! 
EE ee Log|u+y RI. 


F(2n)= 
e 2nQn—2 
sip=2n 


5. Si tma est le segment le plus grand mesuré sur les transects. on calcule l'intégrale entre les 
bornes (h. max) puisque ` 


0, (1—F,(u)=0. pour ` Hz Wenn 
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